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Abstract 

Equilibrium statistics of Hamiltonian systems is correctly described 
by the microcanonical ensemble. Classically this is the manifold of all 
points in the A*'— body phase space with the given total energy. Due to 
Boltzmann's principle, e"^ = tr{5[E — H)), its geometrical size is related 
to the entropy S{E, N,- ■ ■). This definition does not invoke any informa- 
tion theory, no thermodynamic limit, no extensivity, and no homogeneity 
assumption, as are needed in conventional (canonical) thermo-statistics. 
Therefore, it describes the equilibrium statistics of extensive as well of 
non-extensive systems. Due to this fact it is the fundamental definition of 
any classical equilibrium statistics. It can address nuclei and astrophysical 
objects as well. All kind of phase transitions can be distinguished sharply 
and uniquely for even small systems. For transitions in nuclear physics 
the scaling to an hypothetical uncharged nuclear matter with an N/Z— 
ratio like realistic nuclei is not needed. 



1 Introduction 

Classical Thermodynamics and the theory of phase transitions of homogeneous 
and large systems are some of the oldest and best established theories in physics. 
It may look strange to add anything new to it. Let me recapitulate what was 
told to us since > 100 years: 

• Thermodynamics addresses large homogeneous systems at equilibrium (in 
the thermodynamic limit N oo\iM/v=p,homogeneous)- 

• Phase transitions are the positive zeros of the grand-canonical partition 
sum Z{T,^,V) as function of e^^ (Yang-Lee-singularities) . As the parti- 
tion sum for a finite number of particles is always positive, zeros can only 
exist in the thermodynamic limit V^|/5.^ — > oo. 
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• Micro and canonical ensembles are equivalent. ^ 

• Thermodynamics works with intensive variables T, P, fj,. 

• Unique Legendre mapping T E. 

• Heat only flows from hot to cold (Clausius) 

• Second Law only in infinite systems when the Poincarre recurrence time 
becomes infinite (much larger than the age of the universe (Boltzmann)). 

Under these constraint only a tiny part of the real world of equilibrium sys- 
tems can be treated. The ubiquitous non-homogeneous systems: nuclei, clusters, 
polymer, soft matter (biological) systems, but also the largest, astrophysical sys- 
tems are not covered. Even normal systems with short-range coupling at phase 
separations are inhomogeneous and can only be treated within conventional 
homogeneous thermodynamics (e.g. van-der-Waals theory) by bridging the un- 
stable region of negative compressibility by an artificial Maxwell construction. 
Thus even the original goal, for which Thermodynamics was invented some 150 
years ago, the description of steam engines is only artificially solved. There is 
no (grand-)canonical ensemble of phase separated and, consequently, inhomo- 
geneous, configurations. It has a deep reason as I will discuss below: here the 
systems have a negative heat capacity (resp. susceptibility). This, however, is 
impossible in the (grand-)canonical theory where c oc (SE)'^ 

^ How does one prove this?: The general hnk between the microcanonical probabiUty 
gS(E,JV) g^j^^ ^j^g grand-canonical partition function Z{T, ^) is by the Laplace transform: 

Z{T,f^,V) = ^°°^dAre-[-E-A^iv-TS(E,iV)]/T 

de dn e-^I'=-Mn-T»(e,n)]/T (3) 

g const. +lin.+quadr. /gs 




If s{e,n) is concave then there is a single point es,ns with 

J_ _ 9S 
T ~ SB 3 
/J _ dS 

T ~ ~ m s 

In the thermodynamic limit {V oo) the quadratic approximation in equ.j^ becomes exact 
and there is a one to one mapping of the microscopic mechanical e = E/V,n = N/V to the 
intensive T,fi. This on the other hand shows explicitely that the intensive variables 
T, n are ill defined in small systems like nuclei. They, as well as free energy F = E — TS, 
or enthalpy H = E + PV , etc., should not be used. 
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1.1 Example: Why one should not use conventional scal- 
ing theory to fix phase transitions in Nuclear Physics 

Still many authors addressing phase-transitions in nuclear physics believe they 
should link these phenomena by finite size scaling to "nuclear matter" in the 
thermodynamic limit. This is the way how phase transitions can be located in 
conventional canonical statistics of ordinary extensive systems with short-range 
interactions. 

However, infinite nuclear matter with the same N jZ ratio as ordinary nuclei 
does not exist. The Coulomb self-energy of such a system would be infinite. 
Therefore the Coulomb interaction between the fragments of a hot nucleus must 
be switched off in such studies. There is no sense in using conventional scaling 
analysis to relate transition phenomena seen in nuclear collisions to anything 
that does not exist in real life. The Coulomb field has an important effect on the 
fragmentation of real nuclei (see e.g. the article by LeFevre at this conference) 
which one must understand by statistical fragmentation simulations. Thus, in 
the conventional approach there is the only choice, either one describes realistic 
but small nuclei without any connection to phase transitions. Or one has a 
scaling theory to something that does not exist in reality: uncharged nuclear 
matter with fixed neutron to proton ratio. 

If anything like phase-transitions exist in real nuclei then only in the sense of 
a new generalized thermodynamics. In the rest of this paper I will describe such 
a generalization which takes Boltzmann's principle serious and avoids the ther- 
modynamic limit. This opens Thermodynamics to the much larger world of non- 
extensive systems. The most prominent example are of course self-gravitating 
astrophysical systems which I will discuss here also. 

2 Boltzmann's principle 

The Microcanonical ensemble is the ensemble (manifold) of all possible points 
in the 6A'' dimensional phase space at the prescribed sharp energy E: 



Thermodynamics addresses the whole ensemble. It is ruled by the topology of 
the geometrical size W{E, N, ■ ■ •), Boltzmann's principle: 



is the most fundamental definition of the entropy S. Entropy and with it micro- 
canonical thermodynamics has therefore a pure mechanical, geometrical foun- 
dation. No information theoretical formulation is needed. Moreover, in contrast 
to the canonical entropy, S{E,N,..) is everywhere single valued and multiple 



W{E, N, V) = eotrd{E - Hn) 




S=k*lnW 



(4) 
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differentiable. There is no need for extensivity, no need of concavity, no need of 
additivity, and no need of the thermodynamic hmit. This is a great advantage 
of the geometric foundation of equihbrium statistics over the conventional defi- 
nition by the Boltzmann-Gibbs canonical theory. However, addressing entropy 
to finite eventually small systems we will face a new problem with Zermelo's 
objection against the monotonic rise of entropy, the Second Law, when the 
system approaches its equilibrium. This is discussed elsewhere QlEI- A fur- 
ther comment: In contrast to many authors like Schrodinger |3] our ensemble 
is not an ensemble of non-interacting replicas of the considered system which 
may exchange energy. I do not consider the different ways to distribute energy 
over the different replicas. I consider the manifold of the same system at the 
precisely given energy under all possible different distributions of the momenta 
and positions of its constituents (particles) in the 67V-dimensional phase space. 
The result is then the average behaviour when one does not know the precise 
position and momentum of every particle but only the total energy. 



3 Topological properties of S{E, ■ ■ •) 

The topology of the Hessian of s{E, • • •), the determinant of curvatures of s(e, n), 
determines completely all kinds phase transitions. This is evidently so, because 
^S{E)-E/T ^Yie weight of each energy in the canonical partition sum at given 
T, see footnote H Consequently, at phase separation this has at least two 
maxima, the two phases. And in between two maxima there must be a minimum 
where the curvature of S{E) is positive. I.e. the positive curvature detects 
phase separation. This is of course also in the case of several conserved control 
parameters. 



(i(e, n) 



dedn 



Ai > A2 



Ai A2 



(5) 



^eigenvectors 



Of course for a finite system each of these maxima of S{E, • • •) — E/T have a 
non-vanishing width. These are intrinsic fluctuations in each phase. 



3.1 Unambiguous signal of phase transitions in a "Small" 
system 

Nevertheless, the whole zoo of phase-transitions can be sharply seen and distin- 
guished. This is here demonstrated for the Potts-gas model on a two dimensional 
lattic oi finite size of 50 x 50 lattice points, c.f. fig.(^. 
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Figure 1: Global phase diagram or conture plot of the curvature determinant 
(Hessian), eqn. of the 2-dim Potts-3 lattice gas with 50 * 50 lattice points, 
n is the number of particles per lattice point, e is the total energy per lattice 
point. The Hue (-2,1) to (0,0) is the ground-state energy of the lattice-gas 
as function of n. The most right curve is the locus of configurations with 
completely random spin-orientations (maximum entropy). The whole physics 
of the model plays between these two boundaries. At the dark-gray lines the 
Hessian is det = 0,this is the boundary of the region of phase separation (the 
triangle APmB) with a negative Hessian (Ai > 0, A2 < 0). Here, we have 
Pseudo-Riemannian geometry. At the light-gray lines is a minimum of det(e, n) 
in the direction of the largest curvature (vA„„^-Vdct = 0) and det — 0, these 
are lines of second order transition. In the triangle APmC is the pure ordered 
(solid) phase (det > 0, Ai < 0). Above and right of the line CPmB is the pure 
disordered (gas) phase (det > 0, Ai < 0). The crossing P„i of the boundary 
lines is a multi-critical point. It is also the critical end-point of the region of 
phase separation (det < 0, Ai > 0, A2 < 0). The hght-gray region around the 
multi-critical point corresponds to a flat, horizontal region of det(e,n) ~ 
and VAi~ and consequently to a somewhat extended cylindrical region of 
s(e,n), details see ^ C is the analytically known position of the critical 
point (second order transition) which the ordinary q = 'i Potts model (without 
vacancies) wow/rf have in the thermodynamic limit. 
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3.2 Systematic of phase transitions in the micro-canonical 
ensemble without invoking the thermodynamic limit 

• A single stable phase of course with some intrinsic fluctuations (width) by 
Ai < 0. flere s(e,rt) is concave (downwards bending) in both directions. 
Then there is a one to one mapping of the canonical ^s-^the micro-ensemble. 




o r cl e I ' p a r a in e t e r 



• A transition of first order with phase separation and surface tension is 
indicated by Ai(e, n) > 0. s(e, n) has a convex intruder (upwards bending) 
in the direction Vi of the largest curvature > which can be identified 
with the order parameter Three solutions of 

dS 



1 

T 

^ T 

determine the intensive temperature T 
Tv. 



dE , 
^ 'dN 

— 1/(3 and the chemical potential 



ft 
o 

5-1 

d 

0) 




Ql O2 03_ 

OTcier parameter 
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In the thermodynamic hmit the whole region {01,03} is mapped into a 
single point in the canonical ensemble which is consequently non-local in 
o. I.e. if the curvature of S{E, N) is Ai > both ensembles are not 
equivalent even in the limit. 

• A continuous ("second order") transition with vanishing surface tension, 
where two neighboring phases become indistinguishable. This is indicated 
by a line with Ai = and extremum of Ai in the direction of order param- 
eter va=o ■ VAi = 0. These are the catastrophes of the Laplace transform 
E -^T. 

3.3 CURVATURE 

We saw that the curvature (Hessian) of S{E,N, ■ ■ •) controls the phase transi- 
tions. What is the physics behind the curvature? For short-range force it is 
linked to the interphase surface tension. 

Table 1: Parameters of the liquid-gas transition of small sodium clusters 
(MMMC-calculation 0) in comparison with the bulk for a rising number Nq 
of atoms, Nsurf is the average number of surface atoms (estimated here as 

J2 ^cluster) of cluStCrS with Ni > 2 together, a /Ttr = As^ur-/ * Nq/ Nsurf 

corresponds to the surface tension. Its bulk value is adjusted to agree with 
the experimental values of the as parameter which we used in the liquid-drop 
formula for the binding energies of small clusters, c.f. Brechignac et al. ,6j, and 
which are used in this calculation ^ for the individual clusters. 





No 


200 


1000 


3000 


bulk 




Ttr [K] 


940 


990 


1095 


1156 




qiat [eV] 


0.82 


0.91 


0.94 


0.923 


Na 


^boil 


10.1 


10.7 


9.9 


9.267 






0.55 


0.56 


0.44 






Nsurf 


39.94 


98.53 


186.6 


00 




(J /Ttr 


2.75 


5.68 


7.07 


7.41 
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Figure 2: MMMC simulation of the entropy s(e) per atom (e in eV per 
atom) of a system of No = 1000 sodium atoms at an external pressure of 1 atm. 
At the energy e < ei the system is in the pure liquid phase and at e > 63 in 
the pure gas phase, of course with fluctuations. The latent heat per atom is 
qiat = 63 — 61. Attention: the curve s(e) is artifically sheared by subtracting 
a linear function 25 + e * 11.5 in order to make the convex intruder visible. 
s(e) is always a steeply monotonic rising function. We clearly see the global 
concave (downwards bending) nature of s(e) and its convex intruder. Its depth 
is the entropy loss due to additional correlations by the interfaces. It scales 
oc N^^/^. From this one can calculate the surface tension per surface atom 
Usurf/Ttr — Assur/ * No/Ngurf- The doublc tangent (Gibbs construction) is 
the concave hull of s(e). Its derivative gives the Maxwell line in the caloric 
curve T(e) at Ttr- In the thermodynamic limit the intruder would disappear and 
s(e) would approach the double tangent from below. Nevertheless, even there, 
the probability oc 6^" of configurations with phase-separations are suppressed 
by the (infinitesimal small) factor e~^ relative to the pure phases and the 
distribution remains strictly bimodal in the canonical ensemble. The region 
61 < 6 < 63 of phase separation gets lost. 
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3.4 Heat can flow from cold to hot 

Potts — Uodel, q=10 

0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.'^ 

0.008 

"'O 0.006 

+ 

— - 0.004 
I 

CO 
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12.5 
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^-25.0 
-6d.5 
-100.0 

0.2 0.4 0.6 0.8 1.0 1.2 

e = E/N 

Figure 3: Potts model, {q = 10) in the region of phase separation. At e\ the 

system is in the pure ordered phase, at 63 in the pure disordered phase. A httle 
above ei the temperature T = is higher than a httle below 63. Combining 
two parts of the system: one at the energy e\ + 5e and at the temperature Ti , 
the other at the energy 63 — 5e and at the temperature < Ti wih cquihbrize 
with a rise of its entropy, a drop of Ti (cooling) and an energy flow (heat) from 
3^1: i.e.: Heat flows from cold to hot! Clausius formulation of the Second 
Law is violated. Evidently, this is not any pecuharity of gravitating systems! 
This is a generic situation within classical thermodynamics even of systems with 
short-range coupling and has nothing to do with long range interaction. 
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4 Negative heat capacity as signal for a phase 
transition of first order. 



4.1 Nuclear Physics 

A very detailed illustration of the appearance of negative heat capacities is 
given by d.Agostino et al. [^HI- Here I want to remember one of the oldest 
experimental finding of a "back" -bending caloric curve in Nuclear Physics. 
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Figure 4: Experimental excitation energy per nucleon e* versus apparent tem- 
perature Tapp for backward p, d, t and a together with heavy evaporation 
residues out of incomplete fusion of 701 Mev ^^Si-l-^'^'Mo. The dotted curves 
give the Fermi-gas caloric curves for the level-density parameter a = 6 to 12. 
(Chbihi et al. Eur.Phys.J. A 1999) 
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4.2 Atomic clusters 



Here I show the simulation of a typical fragmentation transition of a system 
of 3000 sodium atoms interacting by realistic (many-body) forces. To compare 
with usual macroscopic conditions, the calculations were done at each energy us- 
ing a volume V{E) such that the microcanonical pressure P — §^/§§ — latm. 
The inserts above give the mass distribution at the various points. The label 
"4:1.295" means 1.295 quadrimers on average. This gives a detailed insight into 
what happens with rising excitation energy over the transition region: At the 
beginning (e* ~ 0.442 eV) the liquid sodium drop evaporates 329 single atoms 
and 7.876 dimers and 1.295 quadrimers on average. At energies e >~ leV the 
drop starts to fragment into several small droplets ("intermediate mass frag- 
ments") e.g. at point 3: 2726 monomers, 80 dimers, ^--^5 trimers, ~15 quadrimers 
and a few heavier ones up to 10-mers. The evaporation residue disappears. This 
multifragmentation finishes at point 4. It induces the strong backward swing of 
the caloric curve T{E). Above point 4 one has a gas of free monomers and at 
the beginning a few dimers. This transition scenario has a lot similarity with 
nuclear multifragmentation. It is also shown how the total interphase surface 



area, proportional to N'^f'^^ = X)i ^t'"^ with Ni>2 [Ni the number of atoms in 
the ith cluster) stays roughly constant up to point 3 even though the number 
of fragments (A^/v = X^i) monotonic rising with increasing excitation. 
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Figure 5: Cluster fragmentation 
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4.3 Stars 



Self-gravitation leads to a non-extensive potential energy oc N"^ . No thermody- 
namic limit exists for E /N and no canonical treatment makes sense. At negative 
total energies these systems have a negative heat capacity. This was for a long 
time considered as an absurd situation within canonical statistical mechanics 
with its thermodynamic "limit" . However, within our geometric theory this is 
just a simple example of the pseudo-Riemannian topology of the microcanon- 
ical entropy S{E^ N) provided that high densities with their non-gravitational 
physics, like nuclear hydrogen burning, are excluded. We treated the various 
phases of a self-gravitating cloud of particles as function of the total energy 
and angular momentum, c.f. the quoted paper. Clearly these are the most 
important constraint in astrophysics. 

Retrieve 



Contour plots and density profiles 
of a rotating, self-gravitating 
W-body system showing the 
formation of a stable double cluster 
(left) and an unstable ring (right) at 
different energies. The 
double-cluster structure illustrates 
the spontaneous breaking of 
rotational symmetry at intermediate 
energy and high angular 
momentum. 

Read the Article 
Back 



Figure 6: Phases and Phase-Separation in Rotating, Self-Gravitating Systems, 
Physical Review Letters- July 15, 2002, cover-page, by (Votyakov, Hidmi, De 
Martino, Gross 
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Figure 7: Micro canonical phase-diagram of a cloud of self-gravitating and ro- 
tating system as function of the energy and angular-momentum. Outside the 
dashed boundaries only some singular points were calculated. In the mixed 
phase the largest curvature Ai of S{E,L) is positive. Consequently the heat 
capacity or the correspondent susceptibility is negative. This is of course well 
known in astrophysics. However, the new and important point of our finding 
is that within microcanonical thermodynamics this is a generic property of all 
phase transitions of first order, whether there is a short- or a long-range force 
that organizes the system. 
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5 Conclusion 



Entropy has a simple and elementary definition by the area e^i^''^'"') of the 
microcanonical ensemble in the 6N dim. phase space. Canonical ensembles are 
not equivalent to the micro-ensemble in the most interesting situations: 

1. at phase-separation ( >-heat engines !), one gets inhomgeneities, and a 

negative heat capacity or some other negative susceptibility, 

2. Heat can flow from cold to hot. 

3. phase transitions can be localized sharply and unambiguously in small 
classical or quantum systems, there is no need for finite size scaling to 
identify the transition. 

4. also really large self-gravitating systems can be addressed. 

Entropy rises during the approach to equilibrium, AS* > 0, also for small mixing 
systems, i.e. the Second Law is valid even for small systems 

With this geometric foundation thermo-statistics applies not only to ex- 
tensive systems but also to non-extensive ones which have no thermodynamic 
limit. For the application to Nuclear Physics I believe one should not define 
phase transitions by scaling arguments towards a non-existing "nuclear mat- 
ter" . This may overlook the non-extensivity of realistic nuclear systems. The 
Coulomb field of a fragmenting nucleus has an important influence on the mass 
and charge distribution of multifragmentation. 
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